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Abstract. We discuss the calculation of Casimir forces between a collection of 
A-dielectric spheres. This is done by evaluating directly the force on a sphere 
constructed from a stress tensor, rather than an interaction energy. Two and three 
body forces between the spheres are evaluated for setups of two and three sphere 
systems respectively. An approximate large- A limit is also obtained for the functional 
dependence on the number of spheres. 
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1. Introduction 

In a recent set of papers Pfl 121 El S] , Casimir interaction energies have been evaluated 
for collections of compact objects interacting with different force carrying fields 
(electromagnetic and scalar). The approach taken has been to evaluate a suitable energy 
functional integral using a T-matrix, whereby an interaction energy can be deduced, 
normalised with respect to their energy when separated at infinity. This is a clean and 
appealing way of evaluating Casimir energies once the effective action and integration 
measures are known. See also [3 [6] for related earlier work on spherical geometry issues. 

Looking more generally, one may try to evaluate Casimir interactions (forces 
or energies) with definitions that are not necessarily equivalent. For example, the 
Minkowski stress tensor is not compatible with the Lorentz force law; it is rather the 
standard vacuum expression that is. The fact that the definitions lead to different forces 
means they are amenable to experimental testing; three body forces can be measured 
between dielectric spheres [7] and thus discernible outcomes can be tested for the best 
candidate theory Typically the differences will show up in both the scale of the forces 
and the higher order curvature corrections. 

In this talk I summarise recent work we have done on calculating the Casimir force 
on a single sphere, in an A^-sphere setup [8]. By using a multiple scattering approach 
to evaluate essentially the classical Green's function of the configuration, we are able 
to evaluate the force directly on the sphere. It is given in terms of Mie scattering 
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coefficients together with translation matrices that map TE and TM modes between 
different scattering centres. This is calculated explicitly for the case of two and three 
sphere setups (with overall evaluation being performed in Maple). For the case where 
we have a large number of spheres that interact weakly, we can form a new coupling 
constant for the perturbative expansions and look at the dependence on N rather than 
the details of the configurational setup. 

2. The iV-Sphere configuration 

The problem we are addressing is how to calculate the force on a particular sphere 
as a result of all the interactions with the remaining spheres in a particular static 
configuration. Pictorially, the setup we have for the configuration of iV-spheres is shown 
in Figure [TJ We shall take sphere-1 to be located at the coordinate origin. 




Figure 1. The iV-Sphere system consists of N dielectric spheres of radii 
R[l], ■ ■ ■ , R[N] each centred on N separate coordinate systems Ei, ■ • • ,Ejy, all 
contained in a background dielectric. 



The Casimir force on sphere-1 (in the j-direction) due to the effects of the iV-sphere 
system of differing material properties is given by 

F j (l\N-l)=[ d 3 xViT ij (x). (1) 
Jb 2 

The stress tensor is given by the standard vacuum expression (which is consistent with 
the Lorentz force law [9]) 

Tyix) = Ei{x)Ej(x) + Bi^Bjix) - ^(|E(x)| 2 + |B(x)| 2 ), (2) 

where x G B 2 and it is understood that we are taking the limit for the initial and final 
points. We then need to evaluate the scattering correlation functions (whilst dropping 
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the direct modes of propagation), viz 

poo poo 

limE i (o;)E i (y)= / / dudw'(E™ \x; u)^f (y; a/)), (3) 
^ x Jo Jo 

and similarly for magnetic fields. To construct the scattering two point function we 

write the fields in a mode decomposition [TU] of spherical vector wave functions that are 

centred on each sphere centre. Then by applying the standard continuity equations at 

each of the spheres surfaces, one can calculate the out modes in terms of the in modes 

and scattering (Mie) coefficients. Care must be taken when evaluating effectively the 

classical scattering Green's tensor for the N spheres. It is necessary to choose the centre 

of the sphere where the scattering takes place last as the coordinate origin to evaluate 

the derivative, whilst translating this to the centre of the sphere where the force is being 

evaluated (so as to be able to evaluate the eigenfunctions on the sphere). Assuming that 

the background in which we are evaluating this is filled with quantum noise such that 

the noise-current two point function is non-zero we find for the iV-body force on sphere-1 

(suppressing the 5*0(3) indices) 

F[1\N - 1] = - (-1)*— / ci^coth(^/ir i3 T)(l|[a 1 (w J R[l]) 
x ^^(rM) • a\ujR[7\) J>«(r[i, j}) ■ a> 

i=2 j=2 
N 

x^Vr^^rfj, l])}j(kR[l})h + (kR[l})W(ujR[l])\l).(4) 

Here, a l (ujR{i]) are the Mie scattering coefficients in the £0(3) basis for sphere i with 
radius R[i]. In addition A l ^{r[i,j]) are the translation matrices mapping the modes 
between spheres i and j, and the vectors |1) give the truncation in the L angular 
momentum quantum number (leading to a multipole type expansion). Note the explicit 
form of the translation matrices involve exponentials of the inter-sphere separations [8] 
and thus it is the total path length that plays the key role in understanding the variables 
of the system. 



3. Two and Three Sphere Forces 

For simple setups we can evaluate Equation Q perturbatively. In the case of two 
spheres aligned along the z-axis one can calculate the force at T = as a multipole 
expansion and at T > by residues. An example plot is given in Figure [2} In a similar 
fashion when the configuration consists of three spheres we can evaluate the three-body 
force (where the three two-body forces have been subtracted). In Figure [3] a plot of the 
potential (derived from the integral of force with respect to the appropriate separation 
vectors) on one of the spheres as a function of one of the others position (the other being 
held fixed) is given. 
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Figure 2. A plot of the inter-sphere retarded force between two dielectric spheres 
in the empty vacuum at T = 0°K. Here, the relative dielectric permittivity of the 
polystyrene spheres is t\ — 62 — 2.6 . The angular momentum series representation 
has been truncated at L truncate = 3. 
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Figure 3. A plot of the inter-sphere retarded potential (in units of Hc/Att) between 
three dielectric spheres in silicone fluid with €b = 2.2 at T — 2$S°K. Spheres 1 and 2 
are held fixed along the z axis 10i?[l] apart. Here, the relative dielectric permittivity 
of the polystyrene spheres is t\ = ti = £3 = 2.6 . The angular momentum series 
representation has been truncated at L trU ncate — 1- 

4. A Large N Weak Coupling Limit 

Consider the case of iV identical spheres arranged in some fashion as iV — > 00, 
a 1 = a — >• 0, whilst A := Nas is held fixed (a = asuj 3 R 3 /c 3 ). Also assume that 
there is some representative separation, s, between the spheres and a corresponding 
orientation. In this case there are (N — 1)! irreducibly connected diagrams that 
contribute approximately equally to the iV-body scattering. The expression for the 
associated potential reduces to (considering the case T = to be specific) 

V[l\N-l]n±(-l) N ^(N-l)\ ^ dXe~ x [aA { X,s)] N (5) 

A s Jo 
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where A is a simple polynomial in X obtained from the translation matrices A after 
extracting the exponential prefactor. Taking the limit N — > oo we can extract the 
dominant L = 1 term which reduces to (using Stirling's approximation) 

p -N D37V 

V[l\N-l]~±(-l)»hc^\»£-. (6) 

The use of this formula would be in adding an additional sphere to the system and 
measuring the resulting oscillation of the force compared to the absolute force before the 
sphere is added. Alternatively it could be used for quantifying the error of a truncated 
series representation for large numbers of weakly interacting particles. 



5. Conclusions 



In this talk I have summarised recent work we have done on calculating Casimir forces 
between spheres using a multiple scattering approach. This has been done at both 
zero and finite temperature. The total closed path length being always greater than 
the respective length scales of the associated radii is what enables the perturbative 
evaluation to be an accurate approximation. For large N and weakly scattering spheres, 
we have deduced a scaling type formula for the force based on an average set of properties 
for the configuration. This could be useful for forces where we want to understand the 
behaviour as a function of N rather than the detailed configuration. 
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